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In this talk [1] we report on the recent progresses on IR logarithms resummation for the 
Thrust distribution in e^e~ collisions. Using renormalisation group (RG) evolution in 
Laplace space, the resummation of logarithmically enhanced corrections is performed 
to next-to-next-to-leading logarithmic (NNLL) accuracy. To combine the resummed 
expressions with the fixed-order results, we derive the log(i?)-matching and 7?-matching 
of the NNLL approximation to the fixed-order NNLO distribution. 

1 Introduction 

Event-shapes are observables which measure the geometrical properties of energy-momentum 
^ I flow in a hadronic final state. They have been precisely measured over a broad range in 

energies at electron-positron colliders. The event-shape distributions allow for a detailed 
probe of the dynamics of QCD and especially for a precise determination of the strong 
coupling constant ag. Owing to their infrared and collinear safety, they can be computed 
systematically in perturbation theory. The fixed-order description, based on a power series 
expansion of the distribution in the strong coupling constant, is reliable over most of the 
kinematical range of the event-shape. In the dijet limit, which is attained for the thrust 
variable [2] as T -> 1, the convergence of the fixed-order expansion is spoilt by large loga- 
rithmic terms log(l — T) at each order in the strong coupling constant, thus it necessitates 
a resummed description. During LEP times, precision studies of a standard set of six event- 
shapes were based on the combination of fixed-order NLO calculations [3-10] with NLL 
resummation [11-13]. To avoid the double counting of terms, both expansions need to be 
matched to each other according to matching procedures such as the R and log(i?) schemes 
[14]. In the recent past, substantial progress was made both on the fixed-order and the 
resummed description of event-shapes. Following the development of new methods for cal- 
culations of QCD jet observables at NNLO [15], the NNLO corrections to e+e~ — > 3 jets 
and related event-shape observables were computed [16-21]. More recently, in the context 
of Soft-CoUinear-EfFective-Theory, the resummation for thrust [24, 25] and the heavy jet 
mass [26] beyond NLL has been performed and applied for a precise determination of as, 
and the framework for the resummation of the jet broadening distributions has been out- 
lined [27,28]. In these calculations, the C(q;^) soft corrections were determined only up to a 
constant term by exploiting the renormalisation group invariance of the cross section. Such 
term is also needed to unambiguously match the resummed distribution to the NNLO result 
in the R scheme. In this talk we report on the direct computation of these corrections and 
we provide a new resummed formula. Finally we match the latter to the existing NNLO 
prediction comparing two different matching schemes. 
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2 Fixed-order and resummed distributions 

The differential thrust distribution in perturbation theory is numerically known at NNLO [17, 
20]. At a centre-of-mass energy Q and for a renormalisation seale /i it reads 



(r) + r 2 

,dC 



+ (1) 



where we defined 



"^ = 2^' "-^g' 

and where a is the total perturbative hadronic cross-section for e^e~ — >■ hadrons. The 
explicit dependence on the renormalisation scale is given by 

= — (r) + 2/3o log(x2)_(r), (3) 

^(->^.)= ^(r)+2 1og(x^)(2/3o^(r) 

+ 2A^(r)) + (2/?olog(4))^^(r). (4) 

The QCD /3-function is defined by the renormalisation group equation for the QCD coupling 
constant 

— - = -as{^)\ Po^ —lit + ...]. (5) 

The normalised thrust cross-section is then defined as 

RT{r)^'-l''-^^eir-r')dr', (6) 

where a is the total cross section for e+e^ — > hadrons. In the two-jet region the fixed-order 
thrust distribution is enhanced by large infrared logarithms which spoil the convergence of 
the perturbative series. The convergence can be restored by resumming the logarithms to 
all orders in the coupling constant. The matched cross section can in general be written as 

Rt{t) = Cias)J:iT,as)+DiT,a,), (7) 

where 



C(a,) = l + 5^Cfca^-, (8) 

oo n+1 

logE(r,a,) Gnrnas'^L"' 
71—1 rn— 1 

= LgiiasL) + g2{asL) + '^p^g^i^asL) + (9) 

TT 



LCWSll 



2 



where L = log(l/T). The function gi encodes aU the leading logarithms, the function g2 
resums all next-to-leading logarithms and so on. The constant terms C'i are required to 
achieve a full N^+'LL accuracy. D(t, as) is a remainder function that vanishes order-by- 
order in perturbation theory in the dijet limit t — > 0. 

In view of matching the NNLL resummed distribution to the NNLO fixed order prediction 
using the i?-matching scheme, we need to include the logarithmically subleading terms C2, 
C3 and G31 in the expansions (8), (9). 

The resummation of the thrust distribution beyond NLL was first achieved in [24] using 
an effective-theory approach and revisited in [29] , where the full analytic expressions for the 
0{a^) constant term C2 and the coefficient G31 were derived. The 0{a^) constant term 
C3 is currently unknown, and a numerical estimate is given in [29] together with the full 
analytic expressions of the functions gi{asL). 



2.1 Factorisation and Resummation 

Factorisation properties of event-shapes have been widely studied in the literature [34-36]. 
Referring to Fig. 1 we recast the cross section (6) as 

X J dwS (^^,asin)^ e(gV -P-P- wQ) + 0{t), (10) 

where we neglected terms of order 0{t) which are absorbed in the remainder function 
D(r, Us). We use the integral representation of the 8-function 

e(gV -p-p-wQ)^ — [ — e''^Q'e-'"='e-^^'e-^^"«, (11) 



27ri ./c ^ 



and the Laplace transform to recast Eq. (10) as 



where we set ~ vQ^ and Nq — e"^'^ . The soft subprocess S {Nq/NQ / ^, as{fi)) describes 
the interaction between the two jets of hard collinear particles through soft gluon exchange. 
It can be therefore defined in a gauge invariant way as a correlator of Wilson lines 

^ " ^/ ^^^^"'°'^E^O|^^(0)^»(0)|fce.fc)Jcnt(r.g)(fce.fc|M/„(G)W^fi(0)|0), 

(13) 



where we defined Ts = w/Q. Wn and Wn arc Wilson lines 



Wn{y) = Pexp ^g J dsn ■ A{ns + y) 



(14) 
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Figure 1: Leading regions in dijet factorisation. 



describing the eikonal interaction of soft gluons with the fast moving quarks along the light- 
Hke directions and fi'' respectively. A{ns + y) in eq. (14) denotes the gluon field in QCD. 
The sum runs over the final states |fceifc) involving k soft particles whose phase space is 
constrained according to the thrust trigger function Jcut{TQ^) ■ Both soft and soft-collinear 
contributions are encoded into the soft subprocess. The two-loop expression was computed 
analytically in [29] by performing direct phase-space cuts. The results are in agreement 
with those presented in [32,33]. The collinear subprocess J [J) describes the decay of 
the jet-initiating hard quark (antiquark) into a jet of collinear particles moving along the 
{n'^) direction. It is therefore an inclusive quantity which can be found in many other 
relevant QCD processes such as deep inelastic scattering and heavy quarks decay [22,30,37]. 
The short-distance hard function H {Q/ii,as{n)) = 1^ (Q/Mj Q;s(/x)) [^ takes into account 
the hard virtual corrections to the quark-antiquark production subprocess. It is free of large 
logarithms and it can be generally defined such that Ec^. (10) reproduces the fixed-order 
cross section up to power suppressed terms. 

Using the Renormalisation Group evolution of the soft and collinear subprocesses [29,30,37], 
we can recast eq. (12) as 



exp< — 2 



1 du 



£(1 U 



«2Q2 k 



2 A{as{k^)) + B{asiuQ^)) 



(15) 



where the two coefficients A{as) and B{as) can be computed in perturbation theory. The 
coefficient A{as) reads 



A{as) = rcusp(as) - /3(as) 



da. 



(16) 



where Fcusp (cts) and rsoft(as) are the cusp and the soft anomalous dimensions respectively. 
The former, together with the coefficient B{as) can be extracted from the asymptotic limit 
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of the Pqq{as,z) splitting function [31,38] as 2 — > 1 

P,,(a„ z) ^ 2^p^ + 2B{as)S{l - z) + ... (17) 

l-L ^ Z) + 

The integration countour in eq. (15) runs parallel to the imaginary axis on the right of all 
singularities of the integrand. From eq. (15) we see that the u-integral in the exponent is 
regularised by the lower bound ^ . Such a bound acts as an infrared regulator which prevents 
the strong coupling constant from being evaluated at non-perturbativc scales (< Aqcd)- 
Then, the contour in eq. (15) should be set away from all the singularities (in particular 
from the Landau pole). Nevertheless, for resummation purposes we can set the contour on 
the left of the Landau singularity since it would contribute with a non-logarithmic effect 
suppressed with some negative power of the center-of-mass energy scale. The inversion of 
the Laplace transform can be performed analytically by using the residue theorem as shown 
in [11,29] and results in 

log igi(A) + 32 (A) + ^/3o53(A) + (^)' Ggi log - , 

(18) 
where 

(71(A) = /i(A), 

32(A) = /2(A) - logr(l - /i(A) - A/{(A)), 

53(A) = /3(A)+(^/{(A) + lA/r(A)) (^(°)(1-7(A))2-V^«(1-7(A))) 

+ /^(A)^(°)(l - 7(A)) + Cf/13o {ie (3/2 - 7^) - ^rVe) • (19) 

The functions fi{X) as well as the constants Ci, C2 and G31 are defined in [29], while the 
C3 constant term is still analytically unknown. We fit the latter numerically using the 
fixed order Monte Carlo parton-level generator EERAD3. The fit is performed by subtracting 
the 0(0;^) logarithmic structure from the fixed-order result and taking (numerically) the 
asymptotic limit t 0. 

EERAD3 is run with a technical cutoff yo = 10~^ which affects the thrust distribution below 
To ~ VyO' This forbids us from probing the far infrared region and we perform the fit for 
values of r larger than tq. Numerical fixed order results are obtained with 6 x 10^ points for 
the leading colour contribution and 10^ points for the subleading colour structures. Because 
of the presence of large fluctuations in the Monte Carlo results, each color contribution is 
fitted separately over an interval where the distribution is stable and the different results 
are combined to find the numerical value of C3. As an alternative approach we first sum 
up all the color contributions and then fit C3. We consider the difference between the two 
approaches as a systematic error and as final result we obtain 

C3 -1050 ± 180(stat.) ± 500(syst.) . (20) 

Considering that there is no statistical correlation between different bin errors, as a different 



Rr{r)=(l + Y^C\[^] jexp 
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possible estimate of the systematic un- 
certainty due to the sizeable fluctuation, 
we varied the fit range observing that 
it does not alter the result in any sig- 
nificant way outside the quoted system- 
atic error margins. In Fig. 2 we vary 
the value of C3 within its error band 
and we study its impact on the distribu- 
tion. We observe that the numerical im- 
pact of C3 on the distributions is less than 
1.5%o and it is therefore completely neg- 
ligible compared the other theoretical un- 
certainties such that the large relative er- 
ror range is tolerable for all practical pur- 
poses. 



2.2 Matching to fixed-order and numerical results 

In this section we match the obtained rcsummed distribution (18) to the NNLO fixed order 
prediction. The matching formalism must avoid double counting and allow to access theo- 
retical uncertainties. We compare the i?-matching and log(i?)-matching scheme described 
in [14]. 

In Figure 3 we compare the differential cross section of the new matched NNLL-hNNLO 
results with the old NLL+NNLO derived in [42]. The modification due to the resummation 
is sizable, leading to a 8% increase of the distribution around the peak region. The effect 
of the additional resummed subleading logarithms becomes progressively less important 
towards the multijet region, where the increase is nevertheless of about 5%. It is interesting 
to note that the matching of NNLO with NNLL resummation shifts the pure NNLO result 
also in the multijet region (Figure 4). This was not the case for NLL+NNLO, for which the 
impact of resummation in the region of large r was negligible. This is another sign of the 
importance of the NNLL contribution. 

The renormalisation scale dependence, which was observed to increase from pure NNLO 
to NLL+NNLO [42,45] because of a mismatch in the cancellation of renormalisation scale 
logarithms, is obtained by varying 0.5 < <2. It decreases at NNLL+NNLO by 20% in 
the peak region compared to NLL+NNLO. The magnitude of the scale uncertainty varies 
between 4% in the 3-jet region and 5% around the peak. In Figure 5 wc compare the 
i?-matching and the log(_R)-matching scheme predictions at NNLL+NNLO. The difference 
between the two matching prescriptions is tiny and lies well below the scale uncertainty. This 
implies a very good stability of the theoretical predictions under variation of the matching 
scheme. 

One further source of arbitrariness is the choice of the logarithms to be resummed. In 
fact, it is not clear whether powers of as log(l/T) or powers of e.g. as log(2/T) have to be 
resummed. The origin of this arbitrariness has to do with how much of the non-logarithmic 
part of the fixed-order prediction is exponentiated together with the logarithms. We can 
express this arbitrariness by introducing a new parameter xl , which rescales the logarithms 
as [14]: L ^ L = log (1/(xlt)). 
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Figure 2: Impact of C3 variation 



LCWSll 



6 



1- 


0.5 


"D 




"B 


0.4 


"D 










p 


0.3 








0.2 




0.1 







O 


0.2 


_j 

z 


0.16 


z 
_J 


0.12 


_J 

z 

51 


0.08 




0.04 







;nce 


0.02 


0.016 


■D 

cz 


0.012 


0) 
Q. 


0.008 


Q) 
TJ 


0.004 









0.1 0.2 0.3 

A=(NNLL+NNLO) - (NLL+NNLO) 



0.2 0.3 



























t 








,. ,, I ., , ,. , ,1 .,,.'\ 





0.1 



0.2 



0.3 



0.4 
1-T 



Figure 3: Comparison of the Thrust distributions with NNLL+NNLO and NLL+NNLO 
accuracy. The plot on the top shows the two distributions, with the uncertainty band due 
to scale dependence. The curve in the middle shows the difference between NNLL+NNLO 
and NLL+NNLO normalised to the NLL+NNLO curve. The impact of the resummation at 
NNLL is an increase in the distribution of order 5-8%. The lowest plot shows the absolute 
scale dependence of the two curves. 




Log(R)-matehing NNLL+NNLO 
R-matching NNLL+NNLO 




Figure 4: Comparison of the Thrust dis- 
tribution at NNLO with the matched 
NNLL+NNLO predictions. The contribu- 
tion of NNLL resummation is sizable over 
the full thrust range. 



Figure 5: Comparison of the results ob- 
tained with the i?-matching scheme and the 
log(i?)-matching scheme . The width of the 
curve shows the uncertainty related to the 
scale variation. 



We can estimate the related uncertainty by varying the parameter xl- In Rcf. [14] 
several prescriptions are given on how to set the correct variation range for for different 
observables. For the sake of simplicity and since we are not performing a fit of the strong 
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Figure 6: Dependence on the resummed logarithms, determined by varying the parame- 
ter XL- The left plot shows the change in the xl dependence between NLL+NNLO and 
NNLL+NNLO. The upper plot shows the distributions with the corresponding uncertainty 
band, in the lower plot we compare only the uncertainties. In the right plot the xl depen- 
dence using the two different matching schemes is shown. 

coupling constant, we choose to vary xl within the canonical interval 0.5 < xl < 2, similarly 
to what is chosen to quantify the renormalisation scale uncertainty. This choice is also close 
to the nominal range of variation proposed in [14]. The impact of this variation is shown 
in Figure 6. The left plots show a comparison of the x^-dependence between NLL-I-NNLO 
and NNLL-I-NNLO predictions. The lower plot allows to quantify the reduction of the 
uncertainty due to a variation x^. Apart from the far infrared region, it is observed to 
decrease by 50% in the peak region. The scale-dependence reduction is smaller towards the 
multijet region, where the contribution of the logarithmic part becomes less important. The 
resummation uncertainty at NNLL-I-NNLO varies between 2% and 3%. In the right plot 
the same comparison is made at NNLL-I-NNLO using the i?-matching and log(i?)-matching 
schemes. We observe a similar XL-dependence in both schemes. 

3 Outlook 

The recent results on event shape resummation improve the description of existing experi- 
mental data. In view of future work at high energy linear colliders and precise determinations 
of the strong coupling constant, N^LL predictions for the remaining Event-Shape observables 
are necessary. Moreover, an additional source of uncertainty is due to power-behaved hadro- 
nisation corrections which get large in the dijet region. Currently there is no deep theoretical 
understanding of such corrections which constitute an important source of theoretical er- 
ror. In the past, these were often computed using leading-logarithmic parton shower Monte 
Carlo programs, which turned out to be clearly insufficient [45] in view of the precision now 
attained by the perturbative description. Systematic approaches to hadronization within 
the dispersive model [41,48] or by using the shape function formalism [25,47] are offering a 
more reliable description. Such corrections are quite sizeable at LEP energies (Fig. 7) while 
they are highly suppressed at future linear colliders energies (Fig. 8). In Fig. 7 we show 
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what the power-corrected distribution looks Hke when compared to the pure perturbative 
answer. Non-perturbative corrections are computed with a dispersive model [48] and both 
the mean effective coupling ao and the strong coupling as are obtained by performing a 
simultaneous fit using ALEPH data at Q ~ 91.2 GeV. Such a fit is purely qualitative since 
the correlation matrix is degenerate when only one data set is used. To perform a mean- 
ingful fit, experimental data over a broader range of energies have to be included. We will 
address this issue in a future publication. 




0.1 0.2 0.3 0.4 



1-T 

Figure 7: Theoretical prediction with 
(blue) and without (red) power correc- 
tions compared to ALEPH data. The non- 
perturbative parameter ao as well as the 
strong coupling constant as are fitted to 
experimental data. 
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Figure 8: Comparison of the perturbative 
(red) and power-corrected (blue) distribu- 
tions at a center of mass energy of 500 
GeV. The plot shows how much hadroni- 
sation corrections get suppressed at typical 
future linear collider energies. 
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